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Abstract — We consider the problem of optimal distributed 
beamforming in a sensor network where the sensors observe a 
dynamic parameter in noise and coherently amplify and forward 
their observations to a fusion center (FC). The FC uses a Kalman 
filter to track the parameter using the observations from the 
sensors, and we show how to And the optimal gain and phase of 
the sensor transmissions under both global and individual power 
constraints in order to minimize the mean squared error (MSE) of 
the parameter estimate. For the case of a global power constraint, 
a closed-form solution can be obtained. A numerical optimization 
is required for individual power constraints, but the problem can 
be relaxed to a semidefinite programming problem (SDP), and 
we show how the optimal solution can be constructed from the 
solution to the SDP. Simulation results show that compared with 
equal power transmission, the use of optimized power control 
can significantly reduce the MSE. 

I. Introduction 

In an analog-based distributed sensor network, the sensor 
nodes multiply their noisy observations by a gain and phase 
and transmit the result to a fusion center (FC). The FC then 
uses the sum of the received signals to estimate the parameter 
The key problem in this setting is to design the optimal gain 
and phase multiplier for each sensor in order to obtain the 
most accurate parameter estimate at the FC. Furthermore, 
these multipliers must be updated in situations where the 
parameter is time-varying. Some examples of prior work on 
this type of problem include |T|, fT], fj], [T], \5\. In |T1, 
an orthogonal multiple access channel (MAC) was assumed 
between the sensor nodes and FC. The FC used a best linear 
unbiased estimator to estimate a static parameter and and the 
optimal power allocation with both sum and individual power 
constraints were investigated to minimize the mean square 
error (MSE). A coherent MAC was considered in (2] and a 
linear minimum mean square error estimator was adopted at 
the FC to estimate the Gaussian source. The optimal power 
allocation problem was solved under a total transmit power 
constraint. A phase-only optimization problem was formulated 
in ||3| and the phase of the transmitted signal from different 
sensor nodes was adjusted such that the received signal at 
the FC can be added coherently to optimize the performance 



of a maximum likelihood (ML) estimator In [4J and [5], the 
parameter of interest was modeled as a dynamic process and 
the FC employed a Kalman filter to track the parameter In J4], 
a power optimization problem was formulated to minimize 
the MSE under a sum power constraint and an asymptotic 
expression for the outage probability of the MSE was derived 
for a large number of sensor nodes. Additionally, the problem 
of minimizing MSE outage probability was studied in |51. 

In this paper, we consider a setup similar to 2] and |I5]. 
We assume that the parameter of interest is a dynamic process 
and the sensor nodes coherently amplify (gain and phase) and 
forward their observations of the parameter to the FC. The 
sensors act like a distributed beamformer, but they are also 
forwarding their observation of the background noise along 
with the measured parameter The FC uses a Kalman filter 
to track the dynamic process, and we show how the transmit 
gain and phase of the sensor can be optimized at each time step 
to minimize the MSE of the parameter estimate. We assume 
that the optimized gain and phase is fed back to the sensor 
from the FC at each step, prior to the next measurement. The 
contributions of this paper are as follows: 

1) We find a closed-form solution for the optimal transmit 
gain and phase that minimizes MSE under a sum power 
constraint. This problem was also solved in [4] using the 
KKT conditions derived in (21. However, our approach 
converts the problem to a Rayleigh quotient maximiza- 
tion problem and results in a simpler and more direct 
solution. 

2) The problem of minimizing the MSE under individual 
sensor power constraints is solved by relaxing it to 
a semi-definite programming (SDP) problem, and then 
proving that the optimal solution can be constructed from 
the SDP solution. 

3) For a suboptimal case where the sensor nodes use equal 
power transmission, we derive an exact expression for 
the MSE outage probability. 



II. System Model 

We model the complex-valued dynamic parameter 
first-order Gauss-Markov process: 



as a 



On = aOn-l + Un , 

where the process noise Un has distribution CN{Q,(j'^). As- 
suming the FC and the sensor node are all configured with a 
single antenna, the received signal at the FC is 



Vn 



(1) 



where h„ = [/ii^„, . . . , /ijv^„]^ and hi^n G C is the chan- 
nel coefficient between the ith sensor and the FC, a„ = 
[fli^n, . . . , aN,nY' is the conjugate of the sensor transmit gain 
and phase, H„ = diag{ft,i. „,..., ft,jv,n}, v„ is Gaussian 
measurement noise at the sensors with covariance V = 
E{v„v^} ~ diag {cr^_i, ■ ■ • , cr^.Ar}, and Wn is additive white 
Gaussian noise at the fusion center with variance cr^. The 
channel parameter is defined as 

where „ is complex Gaussian with zero mean and unit 
variance, di denotes the distance between sensor i and the 
FC, and 7 is the path-loss exponent. 

Based on the above dynamic and observation models, the 
standard Kalman Filter is defined by the following quantities: 

. Prediction Step: 9n\n-i = 

. Prediction MSB: P„|„_i = a^Pn-i\n-i + cr^ 

• Kalman Gain: 

P I , h^n 

a-^H„VH^a„ + _P,i|„_]_a^hjih^a„ + 



Measurement Update: 



i|n— 1 



MSB: 



(2) 



III. Minimizing MSB under Global Power 
Constraint 

In this section, we formulate and solve the problem under 
the assumption that the sensor nodes have a sum power 
constraint. The optimization problem is formulated as 



mm 
a„ 



S.t. 



a5Da„ < P„ 



(3) 



max 7 

where D = diagjo-g + cr« n • • • ) ""e + at}, o'e denotes 
the variance of On and Pmax is the maximum sum transmit 
power. From ([2}, to minimize MSB P„|„, we need to maximize 



fc„a^h„ which is calculated as 



P'n I n - 1 ^n ^n^n 



a„H„VH^a„ + P„|„_ia,Yh„h^a„ + 



Thus, the optimization problem Q is equivalent to 

a^h„h^a,j 



max 

Sin 

S.t. a,f Da„ < P, 



a^H„VH«a„ 



(4) 



Denote the optimal solution to problem (HJl as a* . It is easy 
to verify that the sum transmit power constraint should be met 
with equality a*^Da* = Pmax, so that ^ can be rewritten 

as 



max 



a^hr,h^a,i 



aff(H„VH^ + 

s.t. B-n Da„ = Pmax 



-D)a„ 



(5) 



Problem (|5]l maximizes a Rayleigh quotient under a quadratic 
constraint, which results in a simple closed-form solution. If 
we define B = H„VH^ + p "' D, the optimal solution is 
given by 



Pv 



" Y h,^B-iD-iB-ih,,. 
and the optimal value of (|5]l is calculated as 

«1„ iljill„ <1„ 



B 



a;^(H„VH^ 



-D)a* 



h"B^hn 



which is a random variable that depends on the distribution of 
h„. An upper bound for h^B~^h„ is given by 



(a) 



1=1 



(6) 



where (a) follows from B^^ -< (HVH^)-i. Plugging © 
into (|2l), we obtain a lower bound on the MSB: 



> 



1 - 



1 



P I 

n\n—l 



7i\n—l 



1 



P 



n\n~l 



This lower bound can be asymptoticly achieved with Pmax ^ 
00 or crj ^ a1 j, and the corresponding sensor transmit gain 
and phase is 



a„ — 



,) L 



1 



1 



hi^ncrz 1 ■ 



(7) 



From (|7]l, it can be observed that sensors whose product 
\hi^nWy,i is small will be allocated more transmit power. 



IV. Minimizing MSE under Individual Power 
Constraints 

When the sensor nodes have individual power constraints, 
the optimal distributed beamforming problem becomes 



max 
s.t. 



(8) 



where Pmax,i is the maximum transmit power of the ith sensor 
node. Problem dH) is a quadratically constrained ratio of two 
quadratic functions (QCRQ). Using the approach proposed in 
||6|, the QCRQ problem can be relaxed to an SDP problem. 
Introduce a real auxiliary variable t and define a„ = ta„, so 
that problem (jS) is equivalent to 

3^h„h^a„ 



max 

a„.t 



(9) 



s.t. a,j Dia„ < t -Pmax, 

where Dj = diag{0, • • • , 0, (Tg +cr^_,, 0, • 
further rewrite problem (0 as 



« = 1, 



,7V 



max 

S.t. 



, 0}. Then, we can 
(10) 



a,fH„VH^a„+a2<2^1 

n ^t"-n — ^ ^ max, 2; 



1, 



,iv 



Note that the constraints in problem (fTol i already guarantee 
that t 7^ 0, so the constraint t 7^ is removed. 

' h„h» 



Define a„ = [a^,t]^, 



H„VH,f 
0^ 



0^ 











, and Di 



0^ 







problem dTol l can be written in the compact form 



s.t. 



1^ Hrj.arj 



SO that 
(11) 



3-n CIn&n — 1 
Hi 



a^D,a„ < 
Problem ( fTTT i is equivalent to 
tr(AH„) 



1, 



mjix 
A 

s.t. 



(12) 



tr(AC„) = 1 
tr(AD,) < , 
rank(A) = 1 . 



By relaxing the rank-one constraint on A, we convert 
problem (fT2l i to a standard SDP problem: 



max 

A 

s.t. 



tr(AH„) 

tr(AC„) = 1 , 
tr(AD,) < , 
A ^ . 



(13) 



The above problem can be solved in polynomial time using 
the interior point method. Due to the relaxation of the rank- 
constraint on A, the optimal value of problem ( fTSl l provides 
an upper bound for problem (|8). After obtaining the optimal 
solution A*, a rank-one solution a* can be recovered for the 
original problem (|8). In the following, we show that based on 
A* a rank-one solution a* can be constructed such that a* is 
the optimal solution to problem dHJ. 

Defining A;*,^ as the (/,m)th element of A* and A^ as 
the A^th order leading principal submatrix of A* formed by 
deleting the [N + l)st row and column of A*, we propose the 
following theorem: 

Theorem 1. Define the optimal solution to problem ( liiD as 
A*, then A.'^ = aa^ and the optimal solution to problem dSf 

is given by a* = ^ a . 

y ^N+i,N+i 

Proof: We first utilize the strong duality between problem 
( fTsT i and its dual problem to show the property of the optimal 
solution A*. The dual problem of problem (fTjt is given by 
ilJ: 



mm 



s.t. 



(14) 



N 

E 

1=1 



y,T>, + zC„ - H„ > , 
yi,...,yAr,z > . 



It is easy to verify that there exists strictly feasible points for 
problem ( fT3] l and problem (fl4t . For problem (ITji . we can 
construct 

A-'^ — diag{a6, • • • , a6, 6} , 



where < a < mini 



and b 



<''e+<,i E™ 1 a\h,,_,\^al ^+al ' 

For problem (fl4t . we can randomly select ?// > 0, and set 
large enough such that 



> max ■ 



h^h 



2 ' a2 .) 



Then, according to Slater's theorem, strong duality holds 
between the primal problem (fTjt and the dual problem (fl4t 
and we have the following complementary condition: 

tr(A*G*) = , (15) 

where G* = YJ^,^^V*t>^ + z*C„ - H„ and y* and z* 
denote the optimal solution to problem (fT4l i. Due to the special 
structure of D^, C„ and H„, G* can be expressed as 

Gt 



G* = 



'TV 

0^ 





'Ar+l,W+l 



Where G^, = Y.i=i vX^^ + 2*H„VH» - h„h,f and 



G 



Af+l,Ar+i 



Eti yr^'max.z. Since both A* and 



G* are positive semidefinite, Eq. (flSl l is equivalent to 
A*G* = 



Additionally, with consideration of the structure of G*, we 
have 

A^G^ = . 

Define Vg as a set of vectors orthogonal to the row space 
of G^. Then the column vectors of must belong to 
span(VG) and rank(A^) < rank(VG). Given two matrices 
M and N, we have rank(M + N) > |rank(]V[) - rank(N)| 
Is], thus, a lower bound for rank(G^) is calculated as 

rank(G;r) > rank 2/*D, + ^*H„VH,f ^- rank(h„h^) 

= iV- 1 . 

An upper bound for rank(VG) is then given by 



rank(VG) = iV-rank(G;r) 
< 1 . 



(16) 



Since tr(A*H) = h^A^h„ and tr(A*H) > tr(A-'^H) > 0, 
we have 

a;, ^ rank(A^) > 1 . (17) 

Combining Eqs. ( fTSI l and ( fTTI l. we have 

rank(A^) = 1 . 

Define the rank-one decomposition of A^ as A^ — aa^, 
so that the optimal rank-one solution to problem (fTsT i is 



a* — [a-^, Y A^_|_^_^^J^ . 
If the optimal solution of problems dHJ is a* , then 



a,*ff HVH^a* + al 



< tr(A*H) 



where equality can be achieved provided that an optimal rank- 
one solution exists for problem (fTJi . Since tr(A*Di) < 0, 
A* =/: and D.^ >~ 0, then we have A^^^ > 0, otherwise 
tr(A*Di) > 0, which contradicts the constraints in problem 
(l3[ . Based on a*, the optimal solution to problem ^ is given 
by 

1 



and we have 



a„ n„n„ a„ 
a^^HVH^a* + al 



tr(A*H) 



which verifies the optimality of a* . 



V. Equal Power Allocation 

Here we calculate the MSB outage probability of a sub- 
optimal solution in which each sensor transmits with the 
same power. The outage probability derived here can serve 
as an upper bound for the outage performance of the optimal 
algorithm with individual power constraints. The transmit gain 
is given by 



N 



1 



1 



For this suboptimal approach, the MSB is calculated as 

which is a random variable depending on the distribution of 
the channel parameter h„. Define the outage probability as 

Pout = Pr {P„|„ > e}, so that 



Pr 



afH„VH^ae 



</3 



where /? — 

I p 

E = diag 



Pr {af h„h^ae - /3af H„VH^ae < ^a^} 
Pr {h,f (Da^af D - /3E) h„ < pal} , 



D 



— diag ■ 



[hi.n, • • • , hN,n]- 

Define B = Daga^D — /3E, and denote the eigenvalues of 
B as Ai, • • • , Xn, then the random variable h^Bh„ can be 
viewed as the weighted sum of independent chi-square random 
variables AiXi(2). Based on the results in [9], we have 



P. 



N 

E 



u{Xi 



(18) 

where u{-) is the unit step function. If we let ei > • • • > bn 
denote the eigenvalues of E, then from Weyl's inequality ifTOl 
we have the following bounds: 

af D^ae - /3ei < Ai < af D^a^ - l3eN 

-l3eN-i+i < A,; < -l3eN-i+2 , 2 < z < iV , 

-. Since only Ai can be 



where af D^a^ = E,=i NiZ+Jjd^ 
positive, equation < fT8l ) can be simplified as 



Pout 



1- 



1 



Ai >0, 
Ai < . 



Since it is not possible to evaluate the A^ in closed-form, 
the above outage probability expression must be calculated 
numerically. 
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Fig. i; MSE vs. number of sensors for Pmax — 300 or 3000. 

VI. Simulation Results 

To verify the performance of the proposed optimization 
approaches, the results of several simulation examples are 
described here. In the simulation, the distance to the sensors d, 
is uniformly distributed over [2, 8] and the path loss exponent 
7 is set to 1. The variance cTg is set to 1, and the P„|„_i 
is initialized as 1. The MSE is obtained by averaging over 
300 realizations of h„. The observation noise power tr^ ^ is 
uniformly distributed over [0, 0.5] and the power of the additive 
noise at the FC is set to cr^ = 0.5. Two different sum power 
constraints are considered P,nax = 300 or 3000. To fairly 
compare the results under the sum power constraint and the 
individual power constraint, we set Pmax.i = ^-j^- In Fig- E 
the results show that compared with equal power allocation, 
the optimized power allocation significantly reduces the MSE. 
In fact, adding sensors with equal power allocation actually 
increases the MSE, while the MSE always decreases for the 
optimal methods. The extra flexibility of the global power con- 
straint leads to better performance compared with individual 
power constraints, but the difference is not large. The lower 
bound shows the performance that could be achieved with 
Pmax oo. The theoretical and simulated outage probabilities 
of the equal power allocation is presented in Fig |2] The 
results show that the theoretical analysis matches well with 
the simulations. 

VII. Conclusion 

In this paper, we considered optimal distributed beamform- 
ing for an analog sensor network attempting to track a dynamic 
parameter under both global and individual power constraints. 
For the sum power constraint case, we derived a closed- 
form solution for the optimal sensor transmit gain and phase. 
For individual power constraints, we developed a numerical 



Fig. 2: MSE outage probability for equal power allocation vs. 
sum transmit power for TV = 10 sensors. 

optimization procedure that is guaranteed to find the optimal 
sensor gains and phases. We also derived an exact expression 
for the MSE outage probability of a suboptimal scheme in 
which each sensor transmits with equal power Simulations 
were presented to verify the performance of the optimal 
algorithms and the accuracy of the MSE outage probability 
expression. 
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